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Abstract 

Using effective field theory tecliniques we compute the next to leading order Spin(l)Spin(l) 
<^ • terms in the potential of spinning compact objects at third Post-Newtonian (PN) order, including 

o: 

CN . subleading self- induced finite size effects. This result represents the last ingredient to complete the 

Q \ relevant spin potentials to 3PN order from which the equations of motion follow via a canonical 
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formalism. As an example we include the precession equation. 



INTRODUCTION 



There has been significant progress in calculating higher order corrections to the po- 
tentials for spinning compact binaries. The spin potentials to 3PN order, from which the 
0(8182) contribution to the equations of motion (EOM) follow via a canonical procedure, 
were reported in l|] for direct spin-spin interactions, and in [^ from indirect spin-orbit 
effects. These results were presented using the Newton- Wigner (NW) spin supplement arity 
condition (SSC) at the level of the action. The latter procedure was shown to be accurate 

31] using standard power counting techniques. These results were de- 



up to 4PN order in 



s to compute in the 



rived using the NRGR formalism, which appears to provide simple too 

PN expansion |3|, U, |5|, la, [3, la, |9|] , as well as the extremal limit JIOL lllj. More recently we 



computed the spin potentials to 3PN order in the covariant SSC [l3[, and explicitly showed 
the equivalence with our previous results in [l|, l2]- The computations using the NW SSC at 
the level of the action were also re-derived within the NRGR approach in [iJ] , albeit using 



a different choice of metric parameterizations introduced in 15 1. A recent calculation of the 
8182 Hamiltonian to 3PN appeared in [l6|. The results in 16| helped to clarify the neces- 



sity of taking into account spin-orbit effects to compute the 0(8182) contributions in the 
EOM to 3PN when one is working with the SSC at the level of the action. The equivalence 
between ours, and the more traditional approach of 16| , was shown in [2, ll3| . 



In 



13 



we worked within a Routhian formalism originally introduced in [I4I, and devel- 
oped in |3|, Il3|| within NRGR, which incorporates the (covariant) SSC and its conservation 

n n 

upon evolution in a canonical framework to all orders. Within the formalism of [3|, 113[ 
the subleading spin-orbit effect in the EOM is proportional to S^^, the spin tensor in the 
local frame, and contribute at 0(8182) after one reduces spin to a three vector using the 
SSC and takes into account the transformation between the local and global PN frames. 
UM.a.e,. .h^e<,.va,ence between t.ej,dependent .... co.pu.ed i„ Q, Q and Q. 



with those in 



161 ] using ADM techniques 17|, confirms the validity of the new results at 3PN. 



In 13|] we sketched the necessary steps towards computing the next to leading order 
(NLO) 0(8f) contributions to the potential, which is the subject of this paper. We first 
review the Routhian approach, emphasizing the inclusion of higher dimensional operators 
and the preservation of the SSC under time evolution. We break up the calculation into 



three distinct types of contributions: The non-hnear interactions, the finite size effects, and 



the contribution from the extra piece of the Routhian of 13|, which must be included to 
insure that the SSC is preserved upon evolution. As an example we include the precession 
equation to 3PN. In appendix B we show how our result reduces to the geodesic motion 
around a Kerr background in the extreme mass ratio limit. 

II. NRGR AND SPIN EFFECTS IN THE ROUTHIAN APPROACH 

One can introduce a Routhian to describe spin dynamics in a gravitational field as follows 
|3|,U 

^ = - E (^''V^ + \st^ab,U^, + :^R,ea,{x,)SfSf'^ + . . .^ , (l) 



where the ellipses represent non-linear terms in the curvature necessary to account for the 
mismatch between p and u once the SSC is enforced. Since at 3PN order we can consider 
the covariant SSC to be 5*"^^;, = 0, the higher order terms are irrelevant for our purposes. 
There is of course no obstruction to include higher order effects. The overall minus sign is 
chosen to ensure the spinless Feynman rules are not modified and the equations of motion 
(EOM) follow from 

X r ACah J cab 

'nd\ = o, ^ = {s-\n}, ^ = {v,s''''}, (2) 



Sxi" J ' d\ ' ' '' dX 

where the potential is given hj V = —TZ, and 

K,pj = 5^, K,M = c {P^p/^} = o, (3) 

{x^x'^} = 0, {p^/} = ^i^"^.,^'^^ (4) 

{x^', S"*} = 0, K, 5"''} = cu^["5^1^ {P", S''''} = (5) 

with p^ given by V^ = p'^ + ^001^. S"'^, and V^ the canonical momentum. It is easy to show 



the Mathisson-Papapetrou equations 18[ follow from ([2]), and the Riemann dependent term 
in the Routhian guarantees the SSC is preserved upon evolution |l3l ]. 



Notice that our Routhian is similar to the one introduced in 12|, after the replacement 



2m, 



P Qcdoab'^ ^ 



1 Dp'^S'^'u. 



u^ nig dX vii- 



(7) 



The replacement in ([7]) follow from a change of variables. Furthermore, one can show 
that the net effect of this change 12| modifies the spin-gravity coupling as follows ^ 



2<^a^«" 



-^ufSatU'', 



where 



or equivalently (g = 1, 2) 



^ab^gab^:^g4a^b]^ 



Sf = S^v^ 



and 



s^ = s^ + is^- v's^^X - (^' - <s^nvl + .... 



(8) 



(9) 



(10) 



(11) 



From here it easy to show that, when the Routhian is equivalently written in terms of S 
using ([7]), the covariant SSC is conserved, since 



d 
di 



(S'^'u,) = u,{S''', UoiS'^')} + n,^{5»^ S'''}u, + S'^'u, = 0, 






where 



q ^ ^ 

The expression in ( fT2l) follows from the identity 



(12) 



(13) 



«45'^^7^o(5"'')} = o, 



(14) 



due to the fact that the spin algebra in terms of S modifies the expression in ([6]) by 



^ In addition to this change one also generates higher dimensional operators which lead to effects which are 
beyond our interest in this paper. 



shifting rj"-^ -^ t]"-'^ — ^^^^. The inclusion of finite size effects obviously modify the EOM and 
the constancy of the SSC upon evolution is not guaranteed. However, from (TT2|) and (TT4|l 
it is clear that higher dimensional operators describing the internal structure of the bodies 
that are written in terms of S""^ will preserve the consistency of the covariant SSC as shown 
in y. 

In the NRGR formalism, physically relevant higher dimensional operators are those which 
are written in terms of the electric, Eab, and magnetic, Bab, components of the Weyl tensor 
|4, la, lZI, y] . Terms which are proportional to the Ricci tensor, or scalar, can be removed by a 
field redefinition since they vanish on-shell ^. As we mentioned earlier, in order to preserve 
the SSC constraint upon evolution one needs to use S"''' for the higher dimensional terms in 
the wordline action. However, once these terms are written as a functions of the electric and 
magnetic part of the Weyl tensor, it is easy to show that using S"'^ still preserves the SSC, 
as a consequence of the SSC itself and the orthogonality relation, EabU^ = BabU^ = 0. For 
example, the first higher dimensional operator we encounter is the self-induced quadrupole- 
like term, which written in terms of 8°"^ takes the form [l|, |7|, [l3| (g = 1, 2) 

If we now take Les^ and expand 5'^* in terms of 5'°'' using (Q it is easy to show that the 
difference between (fTSll and 

Les^ = ^-^S^S* (16) 



2mmp yy2 

is proportional to {S"'^UbY and therefore can be set to zero, given that it does not affect the 
EOM since it produces a correction which is proportional to the SSC itself. Therefore, we 
have a choice: we can use either flTSj) or flT6|) . as they lead to the same result. In what follows 
we will use (IT6ll . although we will also provide the equivalent result using ( IT5ll in appendix A. 



The operator in flT6l) reproduces the well known LO spin quadrupole contribution to the 
potential 

yfpN = -C^ES^Tph (Si ■ Si - 3Si ■ nSi ■ n) + 1 ^ 2. (17) 

In the case of a rotating black hole Ces^ = l? and this term represents the non-vanishing 
quadrupole moment of the Kerr solution. The coefficient for other compact objects can 



be calculated via a matching procedure. In this paper we will compute the corrections 
from flT6|) to the gravitational potential at 3PN. As we will show below there are no other 
contributions at 3PN from other higher dimensional operators. 

By using the EFT power counting rules it is easy to organize the perturbative expansion 
in a systematic way. To obtain Post-Newtonian corrections one calculates TZ, or the effective 



potential, perturbatively, without imposing the SSC (see [13l] for details). The advantage 
of this approach is that one does not have to worry about complicated algebraic structures. 
The price to pay is the need of a spin tensor rather than a three vector, though once we 
find the EOM via ([2]) we may write our results in terms of three vector and the coordinate 
velocity. For the former a precession equation can be obtained. 

III. NEXT TO LEADING ORDER S'^ POTENTIALS 

The algorithm to calculate potentials in the EFT approach is quite simple. First of all 
we take all of the terms in the Routhian and collects them according to their order in the 
power counting (see [4, LZl] for details of the power counting). Then we draw all possible 
Feynman diagrams at the order of interest. Each diagram is written in terms of a set of 
scalar integrals, and the diagram adds a term to the effective action given by —i J dtV, 



where V is the contribution of that diagram to the effective potential (see 13 1). Throughout 



201 ] for details on 



this section we will suppress the factors of "/ dt" in the diagrams. See 
EFTs. 

A. Feynman rules: spin graviton vertex 

In the weak gravity limit the relevant spin couplings are [ll, uX Il3l | 

rNRGR _ ^ u oik (->o\ 

^IPN — ^ J^iO,k^ , {i-O) 

L^S = ^ {H,,,kS''u= + i^oo,fc^°') , (19) 

L^^S"" = ^ {Hoj,kS"'u^ + H,o,oS'') + -^S'^ {H}Hox, - H^Hoi,k) • (20) 



We refer the reader to [5| for all the other Feynman rules not involving spin. 

B. S^ terms from non linear gravitational effects 

The Feynman rules in fll8lll9ll2Up will contribute to the potentials terms which do not 
arise from either finite size effects or the Routhian term of the form RSS. These non-linear 
gravitational terms contribute from the diagrams shown in ([1]). 



» n 



FIG. 1: Non-linear S'^ terms. The blob represents a spin insertion. 

It is easy to show that the "seagull" diagram (the one not containing the three graviton 
interaction) contributes at 4PN, since (iJoo-f^oi) = 0. (Where bracketed polynomials corre- 
spond to time ordered products and reduce down to retarded Greens function for quadrat- 
ics.) Therefore, for the sub leading non-linear 5*^ effects the only non- vanishing contribution 
comes from the three graviton vertex which resides in the three point function, (HooHoiHoj) . 
This contraction can be easily handled by a short Mathematica routine ^ which can be found 



at 211]. The result from the non-linear effects at 3PN is 



Fig.[I]= ^'^Sl'Si\6^^ - 2n^n^) = z^(Si ■ n)^ (21) 

from which we obtain the potential 



„2 



m2G 



v:; = -^i^{s^-ny. (22) 



jo4 



^ In jlSJ I it has been shown that this diagram can be ehniinated by a different choice of the metric variables. 



C. 5^ terms from the Routhian 



Recall that the Routhain includes a terms quadratic in the spin given by 



LrS'^ — —7. Rdeab(,Xq)S^ S[ 



cdqab^W;^ 



(23) 



whose LO contribution shows up at 2.5PN, 



t2.5PN 



2m 2m 



(24) 



To get a 3PN net effect we need to calculate a diagram similar to Fig. [3b, with the box 
representing now an insertion of ^^52 , to produce a contribution 



{Lll^'^im.Ho.vi)) = z^^^ [3n ■ (v^ x S,)n^ - (v, x S^)'] S^ 



Vr. 



miv^ 



(25) 



Also at 3PN we have a contribution from 



^RS^ = 7 — y^HiifijS"^^ + 2Hoo,ijS^ — S'^^v {Hij^ki + Hki,ij — Hu^kj — Hkj,ii)) S'^ Vc, (26) 
which contracts with a LO mass insertion to account for a 3PN contribution given by, 






(v2 - 3vi) X Si)' + 6n'(vi x Si) ■ n + 3(n x Si)'(n ■ v) 



Sfv,. 
(27) 



Adding both pieces together we end up with the contibution to the potential. 



^3lS=Kl))'5i%e + l-2, 



(28) 



with sl^Id the Si piece of the acceleration in the local frame, 



m2G 



N 



h{i) 



m^iT-^ 



[-3v X Si + 6n(v X Si) ■ n + 3n ■ v(n x Si)] 



(29) 



Notice that we substituted the SSC for the S{ term inside the bracket in fl26|) since it 



multiplies the SSC itself. From here we can use (^ to obtain the EOM from which we get 

^ = ... + (a?^,)XSOxvi + ... (30) 

A similar expression straightforwardly follows by using the acceleration dependent form 
of the extra piece in the Routhian after the redefinition in ([71). 

D. Subleading 5^ terms from finite size eff"ects 

Let us now consider the finite size corrections. Using the power counting [Sj it is straight- 
forward to show that the LO contribution from ( IT6l) scales Les^ ~ yLv^, for maximally 
rotating compact objects, with Ces^ ~ 1. As we mentioned before, this term generates a 
gravitational potential of the form of flTTj) at 2PN. As it is well known, rotating BHs or 
NSs have a quadrupole moments given by Qhh = —aS'^/m {G = c = 1), with m, 5, the 



mass and spin respectively ! 



For a BH we have a = 1, and for NS a ranges between 4 



and 8 depending on the equation of state of the neutron star matter [23|]. Matching this 
result with the effective theory we find Ces'^ = ct; which is consistent with what we expect 
from naturalness arguments. It also tells us that this contribution is enhanced for NSs with 
respect to BHs, since the latter seems to provide a lower bound for Ces^ iii (H6I) ^. 

In order to calculate the 3PN NLO correction due to the finite size of the objects we need 
to take care of two things: first we need to expand (1161) in powers of the relative velocity 
V, and compute all possible S^-orbit diagrams such that the net scaling goes like v^; and 
second we have to make sure we are not missing any other higher dimensional operator 
which could contribute at 3PN order. 

1. All Possible Higher dimensional operators 

Let us start with the higher dimensional operators. At NLO we have a few operators 
that could contribute. Let us start using 5°'' to begin with since its use guarantees the 
preservation of the SSC upon evolution. Furthermore, let us put aside reparameterization 



3 



Possibly such a bound could follow from studying graviton scattering off the finite sized object using the 



dispersion relation techniques developed in [24 1 



invariance (RPI) for the time being . If we define 5" = e'^^'^'^'SbcUd we can write down the 
following new terms in the action 

DaB.dS'^'S'', (31) 

DaE^dS^'S". (32) 

where Bab represents the magnetic part of the Weyl tensor. 

These terms are self-induced effects which can be generated by diagrams in the one-point 
function with two spin insertions. These operators scale as \/IjV° and ^/Lv'^ respectively. 
Both will contribute beyond 3PN since the magnetic component of the Weyl tensor does not 
couple to a LO mass insertion, e.g. {BabHoo) = 0. We may wonder about operators with 
only two derivatives since we could also have 

EbcS'^^UaS'^'ud (33) 

B^bS'^'^S'^Ua. (34) 

Notice however that these operator cannot contribute since S°'''ub = 0, identically. We 
could, however, have chosen to write these operators in terms of S"''', in which case the first 
one does not contribute, since it is proportional to {S°'^UbY and thus always has a vanishing 
contribution to the EOM. On the other hand, the second term in (IM|) would be equivalent 
to 

RabcdS''S'''u''Ue, (35) 

and we immediately recognize this is our extra term in the Routhian of ([1]). Notice that, 
once the SSC is enforced, this term does not contribute to the n-point function. This makes 
fixing its coefficient, by matching to the full theory, ambiguos. Which is to say that the 
coefficient is independent of the underlying theory and must be fixed algebraically. This is 
consistent with the fact that Wilson coefficient is fixed by the consistency of the SSC, in 
this case the covariant choice, as discussed in the previous section. 
Finally we should consider the term 

Lbs'^ = CBs^BcdS"" Sb, (36) 



10 



similar to the electric quadrupole. The contribution of this operator to the potential stems 
from the coupling to an 0{v) mass insertion of the form HQiv\ The LO potential at 2.5PN 
reads 

V^S2 = Cb52— f (n X V2) ■ Si(n ■ Si) + 1 ^ 2. (37) 

It turns out however that this term has a vanishing Wilson coefficient, e.g. Cbs^ = 0, 
due to parity conservation. For instance we can easily show that it vanishes for the case 
of a Kerr BH by comparison with the multipole expansion of the Kerr metric in harmonic 
coordinates^. 



2. NLO corrections to the (spir? )quadrupole-monopole interaction induced by L 



We are then left to compute the subleading corrections due to ( fT6l) . Using the power 
counting rules of NRGR for spinning bodies [Tl] we obtain the following new vertices 

t2PN ^ES'^ tt Qiknjk /QQ^ 



j^3PN ^ C^ES^ 

^^' 2mmp 



^Hoo.jS'^S^^ + S^'^S^'^ (Ha,ojv' - ^Hir^.v^v' + Hu.rv'v^ (40) 
^2 1 \ . 

+ Hoo^i{Hij^i — Hiij — Hji^i) + Hoi^iHoj^i — Hio^iHjQ^i 
+ HooHqo^ij — 2HiiHoojj\ 

where a' = ^ is the particle's acceleration^, and we have discarded terms which are total 
derivatives at 3PN. At this order we must also include diagrams with double graviton ex- 
change. The quadratic term in (l4T!l arises from expanding '"VI '' to second order in the 



4 



There is no S'^ term in the expansion of the hot component of the Kerr metric 25 1 



^ This term follows from integrating by parts the piece proportional to — ^^^ S^^ S^^hij^ai3V°'v^ in the 
expansion of the electric component of the Weyl tensor in (|16|) . This term was incorrectly discarded in 
our previous code as a total derivative. We thank A. Ross for pointing this out to us. 



11 



metric and tetrad perturbation. Notice that the LO finite size effects in the spin^-spin sector 
contributions, shown in Fig. [2], start at 3.5PN and therefore can be ignored. 



O 



-O 



V 

a) 



V 

b) 



FIG. 2: Leading order finite size spin^-spin contributions. The black square represents an insertion 
of the finite size operator 



V 

a) 

.,6 



V 

b) 



V 

c) 



d) 



FIG. 3: Diagrams contributing to the 3PN which do not involve non-linearities. The cross corre- 
sponds to a propagator correction, i.e. a retardation eff'ect. 

Fig. [3] shows 3PN contribution from one graviton exchange, with the last diagram cor- 
responding to the first retardation effect, while Fig. H] gives the contributions form non- 
linearities. Calculating these diagrams is straightforward. The results for each one of the 



12 



diagrams are 



.(1) 
Ami r-^ 



i C^g2 37 77-2 Gn ^^2QikQJk, 

Ami 



Fig. mx = T ^ wlS'^S{\Zn'n^ - 5'^) (41) 



Fig. [3b = _!%!!!^^ V2 • vi5f 5f (3n^n^- - (5^^), 



7?7i r'^ 



Fig. [Sb = _!^i5£!!!^^ fe°^f (3n*n^' - 5'^) - 25f ^f (t;! - 3vi ■ nn^) 
r?7i zr'^ L 

+ 2S\^St ( (3n • wm^vi - v\v{) - Kl{?,n'n^ - 5'^) + (7;i7;^2 " 3v2 ■ nn^v\) - U'^s^ ■ r J 

for the instantaneous one-graviton exchanges, and 

Fig[3ll = -iC^^l^^^^Sl'Sl' [vi ■ V2(<5'^' - ?,n'n^) - ?,{v{n' + i^lri^Vs ■ n 

-3vi ■ n(7;^?7* + v^n^) + v\v^2 + ""Wi - 3vi ■ nv2 ■ n{5^^ - bn^n^)] (42) 



for the correction stemming from expanding the propagator. 

For the non-hnear terms, since the LO finite size insertion is proprtional to -f^oo.jj? we 
need the three point function {HoqHqqHqo) , which we can get from the results for the spinless 



B 



case in [^, and just correct for the worldhne insertions. We end up with 



FigSk = -^d^l9k!i^4n,n, - 5.,)^f ^f , (43) 

Figib = -zCjJ),^|!^(3n,n, - 5ij)Sl'Si\ (44) 

Figit = ^C^El^j^i^n^n, - 6,,)Sl'St, (45) 

FigiH = -z4^),^^(2n.n, - 5.,)5f 5^ . (46) 



As in [13|, we spht the potential at 3PN into two pieces 



Vgo = 4s^^^ [Si'Sf{3n^n^ - ^ - 2Sf ((v^ x Si)'^ - 3(n ■ vi)(n x S^)')] 

+1^2 (47) 



13 



» « 



» * 



a) 



4 

V V 



b) 



V V 

c) 



d) 



FIG. 4: Non-linear finite size contributions to the 3PN spin-orbit potential. 



and 






C 



(1) <^JV^2 



S? 6(n ■ vi 



15 13 3 2 7 2. 

— n ■ vin ■ V2 + — vi ■ V2 - -v^ - ij^i - 2ai ■ r 



/g 21 15 \ 

+ (Si ■ n)^ ( 2 '^'^i ^ ^^2) - yVi ■ V2 - — n ■ vin ■ vg j + 2vi ■ SiVi ■ Si 

— 3vi ■ S1V2 ■ Si — 6n ■ vin ■ SiVi ■ Si + 9n ■ V2n ■ SiVi ■ Si + 3n ■ vin ■ S1V2 ■ Si] 



-C 



(1) "laG 



N 



ES'^ 



2r4 



(l + ^)(S?-3(Si.nn + 1^2. 



(48) 
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E. The 5^ potential in the covariant SSC 

Collecting all the pieces together we end up with the following expression for the spin^ 
potential to 3PN in the covariant SSC 



V, 



3PN 



v,l,^^ + v^o. + vi + v:, + V, 



'3PN 



'SO' 



nl 



2PN 



(49) 



C 



(1) GN'm2 



(1) G]\im2 



[^f ^f (3nV - 6'^) - 25f ((vi X Si)^ - 3(n ■ vi)(n x Si)'^)] 



+C 



ES^ 



2mir^ 



-(Si-n)^ 



9 



8^(6(11 -vi)^ 
o. 21 



15 



2 

Vi ■ V2 



n ■ Vin ■ V2 

15 



13 



7 



-Vi ■ V2 



2 2 2 2 ^ 

n ■ vin ■ V2 ) + 2vi ■ SiVi ■ Si 



2air 



2V 1 ^^ 2 2 

— 3vi ■ S1V2 ■ Si — 6n ■ vin ■ SiVi ■ Si + 9n • V2n ■ SiVi ■ Si + 3n ■ vin ■ S1V2 ■ Si] 

-CS.i0 (S; - 3(S, . n)^) + C%!^ (1 + i^ 

G?^m2 



(S?-3(Si-n)^) 



[Si ■ nf + (a^^i))' 5f + vi X Si • a^^^) + 1^2. 



In appendix B we provide a cross check for this potential by taking the extreme mass 
ratio limit and showing that we reproduce the motion of a test particle in a Kerr background 
as expected. 

F. Divergences and Regular izat ion 



In computing the Feynman diagrams leading to the expression in fH9|) we encounter 
divergences of many sorts. First of all from Wick contractions such as the ones represented 
in Fig. [51 




FIG. 5: Divergent diagrams which renormalize the mass and quadrupole couphngs. 



(g) 



These diagrams however simply renormalize the couplings of our theory, mg and C^g2 



coefficients. The divergences are power like and are thus "pure counter-terms" (see 20| for 



15 



a discussion for non-experts). Thus, no renormalization ffroup (RG) flow is present in tliis 



case, contrary to the case of logarithmic divergences [4, |7|, ll3|. In the case of dimensional 
regularization, these diagrams are automatically set to zero since they involve scaleless 
integrals. Other divergences are found for instance in diagrams such as Fig.H^b, or similarly 
in Fig. [H These divergences occur when a factor of momentum squared in the numerator 
from the three graviton vertex cancels the intermediate propagator line connecting with the 
external worldline. Effectively the resulting diagrams appears as in FigO These represent 
a (pure counterterm) renormalization of the mass and the finite size coefficient Ces^- 

At higher order, tidally induced logarithmic divergences will appear. As shown in |7|, 
the first tidally induced finite size effects starts out at 5PN for maximally rotating compact 
objects, from a higher dimensional operators of the form 

CD2Es^{ft)D^E,kS:S''. (50) 

In the expression of (150|) Cu2Es'2{fJ') is a Wilson coefficients which obeys a RG-type of 

equation, 

d ^ , , m . . 

l^-r^D^ES^W 7- (51) 



dix m. 



Notice that (ISUj) does not contribute to the one-point function, as it is expected from 
Birkoff's theorem. That is due to the fact that D'^h^j^ ~ 0{h'^) on shell, and the linear 
piece in (1501) can be traded, by a field redefinition, for a term proportional to h^ . Therefore, 
(150|) does contribute to the (n+l)-point function and that is how we can match with the 
full theory using scattering amplitudes. This result generalizes the effacement of internal 
structure up to 5PN order^, modulo self-induced effects computed in this paper, to the case 
of spinning bodies. 

IV. THE S^ CONTRIBUTION TO THE SPIN DYNAMICS TO 3PN IN THE CO- 
VARIANT SSC 



From (1491) we can obtain the spin dynamics using ([2]). The result reads 



^ As it was shown in [^ finite size effects enter at 5PN for spinless bodies via terms proportional to E^, B^ 
in the worldline. These terms do not contribute to the one-point function as well and matching can be 
also achieved via comparison of scattering amplitudes in the EFT and full theory sides. 

16 






cu.T'^ + 00^/"") X Si + (a^?l^ X Si) X vi + (vi X Si) x Ai + Bi x Si (52) 



where we included the term from the Routhian in (l30D. and 



,(1) 'm2GN .Q X 
'2PN - ^^s52^^^^n(^J3i ■ nj. 



^OPAT — oLj t 



(53) 



Also 



Ai = vi X uj. 



2PAf5 



(54) 



and 



Bi = %^ [2vi(Si ■ vi) - 6vi(n ■ vi)(Si ■ n) 



(55) 



c2 _ (1) GNm2 

^3PN — ^ES^ 



2mir-^ 



9 21 15 

2(Si ■ n)n ( -(v^ + v^) - — Vi ■ Vg - — n ■ Vin ■ V2 ) + 4vi(vi ■ Sj 



3vi(v2 ■ Si) - 3v2(vi ■ Si) - 6n(n ■ vi)(vi ■ Si) - 6vi(n ■ vi)(n ■ Si) 



9n(n ■ V2)(vi ■ Si) + 9vi(n ■ V2)(n ■ Si) + 3n(n ■ Vi)(v2 ■ Si) 



+ 3v.(n.Vi)(n.Si)]-Cg.^|l 



—^ 3n Si ■ n V^n Si ■ n . 56 

mi / r^ 



A. The precession equation 



The spin EOM in ( |52l) can be transformed into a 



transformation to the NW SSC which we applied in 13 | 



precession form by performing the 



Sr = ( 1 - \^i) Si + ^vi(vi ■ Si) + || [(n X S2) X V2] X Si 



(57) 



with V the velocity in the local frame given by 



GN'm2\ Gn c< , 
vi = I 1 H vi H — n X S2 + 

1" / ij"-^ 



(58) 



and 



Xq > Xq 



2mn 



(vg X Sq) + . . . , 



(59) 
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where the elhpses represent higher order corrections. 

The EOM in terms of S""' takes a precession form. This form was aheady shown for the 

|-Sr = ^7 X Sr (60) 

with (we suppressed the nw label in the spin vector for simplicity) 

uj = uj^ + ul^"" + Bi + ^ (ai(i) + 0. J^ X vi) X vi + 5ut + 5^,^^ (61) 

where ai(i) and Bi are given by ( l29l) and ( l55l) respectively. We also have the 0{Si) pieces of 



the shift in the spin-orbit frequency after ( 1571) 
into the NW spin-orbit precession as shown in 
inM), 



which effectively takes the covariant result 
and hence from (!59|) (see also eq. (74) 



^^fo = ^^^ |n X ( -vi - 6V2J (n X vi) ■ Si + viSi ■ ( -vi - 2v2 

3G>^ (62) 

Also from the inverse of (1571) in the LO finite size term 



2PN _ ^(1) 3GArm2 



^-^l'"" = ^ES^^^ [^(^ ■ Si)v? - n(n . vi)(vi ■ S^)] . (63) 
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V. THE SPIN POTENTIAL TO 3PN IN THE COVARIANT SSC 

For completeness here we re-write the spin-spin and spin^ potentials, including also the 
LO spin-orbit term, 

Gn 






{5'' - 3nV^) ( 5f 5f + ^vi ■ V25f ^f + v'^'v^^Si^Sr - v^v^Si^St 



+ 5f 5f'(^2 - ^'t) + SfSi\v\ - vf)\ + -5f ^2^^' (3vi ■ nv2 ■ n{5'^ - hn'n^) 
+ 3vi • n{v{n'' + Wg'^-') + 3v2 • n{v{n^ + t;^n'') — v\v{ — v\v{) 
+ (3r2'v2 ■ n - v[)Sl^Sf + (3r2'vi ■ n - v[)SfSl^] 

+ [CeI^^ [Si'SfiSn^n^ - ^ - 2Sf ((v^ x S^^' - 3(n ■ vi)(n x S^)')] 



^(1) GNm2 



Sj { 6(n ■ vi)^ - — n ■ Vin ■ V2 + — Vi ■ V2 - -v^ - -v^ - 2ai ■ r 



+ (Si ■ n)^ ( -(v^ + v^) - — vi ■ V2 - — n ■ Vin ■ Vs j + 2vi ■ SiVi ■ Si 

— 3vi ■ S1V2 ■ Si — 6n ■ vin ■ SiVi ■ Si + 9n ■ V2n ■ SiVi ■ Si + 3n ■ vin ■ S1V2 ■ S 

-^'eI'^ (S? - 3(S, . nr) + C« !!|^ (1 + *-;^) (S; - 3(S. • n)=) 



which allows us to obtain the (9(SiS2) (computed in [13|) and 0{S^,2)) contributions 
to the EOM to 3PN in the covariant SSC. The latter is imposed after the EOM are 
obtained via ([2]). Note that, to the order we are working in this paper, we can replace 
the acceleration-dependent term in our Lagrangian by the LO equations of motion, or 
equivalently integrate it by parts. 



Missing in (1641) is the NLO spin-orbit potential. The latter was obtained in J26|, |27|. 
However, recall that our computation is in terms of the spin in the local frame, S""^, whereas 
in 



26 



27| the spin dynamics was calculated in terms of the spin tensor in the PN frame. 



e.g. S^'^. The result can be translated to the local frame by using the tetrad field e" |l3 |. 
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The potential in (1641) therefore completes the computation of spin effects to 3PN. 



VI. CONCLUSIONS 

In this paper we have completed the ca 
3PN. Together with our previous results in 



cu 



ation of the potential quadratic in spins to 
2|, [l3|, the results in this paper complete the 



computation of the relevant potentials to obtain the spin effects in the EOM to 3PN order. 
It is important to note that within NRGR there is no obstacle to go beyond this order using 



the formalism developed in 4, [iSj. In particular, the NNLO corrections should follow with 
relative ease. The radiation at 3PN from spinning binaries remains to be calculated. The 
framework for such a calculation within NRGR was set up in ^. We report on these effects 



m 



28|. 
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APPENDIX A: USING S INSTEAD OF S 

In this appendix we compute the subleading finite size corrections using flTSj) . which we 
will explicitly show give rise to identical results to those found in the body of the paper. The 
calculation is identical, the only exception is that now we use S"'^ in the diagrams. As we 
did in 13|, we can sum up all the diagrams and hence split the finite size potential at 3PN 
in two pieces, as we did before in (j48j) and (147|) . Hence we have again one term dependent 
on 5*, 

VSl = ^5^1^ [5f 5f (3nV - 6'^) - 25f ((vi x Si)^ - 3(n ■ vi)(n x Si)^)] +1^2 

= ^5^1^ [-3(vi X Si)^ + 3 (n ■ (vi X Si))' + 6(n • vi)(vi x Si) • (n x Si)] 
+1-^2 (Al) 
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where we used f lTUl) . and as before at 3PN, 



Vi 



(1) GNm2 
^•^^ 2mir3 
9 



Sj { 6(n ■ vi)^ - — n ■ vin ■ V2 - yVi ■ V2 - -v^ - -v^ - 2ai ■ r 



+ (Sl-n; J -vv^ 



2 I ■-' t 2 I 



vi ■ V2 - — n ■ Vin ■ V2 ) + 2vi ■ SiVi ■ Si 



3vi ■ S1V2 ■ Si — 6n ■ vin ■ SiVi ■ Si + 9n ■ V2n ■ SiVi ■ Si + 3n ■ vin ■ S1V2 ■ Si 
(1) rn2G'iT ( . 4m 



-&^J^^:^^ I 1 



'ES'^ 



2r4 



■m\ 



(S?-3(Si.nf) 



(A2) 



We are still one term short to complete the potential to 3PN. We still need to include 
the subleading corrections due to ([9]) in the LO finite size potential of ([1] 



V^l,, ^ &llJ^^2^SfS{\8'^ - 3n^n^) + 1-^2. 



'2PN 



'^^'^ 2mir^ 



(A3) 



Equivalently we have 



V, 



2PN 



(1) m2GN 



{Si ■ Si - 3Si ■ n5i ■ n) + 1 ^ 2, 



(A4) 



with Si^ = e'^^S^, and 

{^;,5j} = -6^^'(^; + (S,-v>J). (A5) 

The final result for the spin^ potential in the covariant SSC turns out to be the sum of 



each previous computations in 



, (HZD, (HHD and dMl), resulting in 



^ipN — ^3PAf + Ks« + ^S + Kl + ^2Pn{^) 



(A6) 



C%^^ [-3(vi X Si)2 + 3 (n ■ (vi X Si))' + 6(n ■ vi)(vi x Si) ■ (n x Si)] 



^ ^(1) GNm2 
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13 



S'i ( 6(n ■ vi)^ - yii ■ vin ■ V2 + yVi ■ V2 - -v^ - -v^ - 2ai ■ r 



+ (Si-n)^(^(vi+v^ 



vi ■ V2 - — n ■ Vin ■ V2 ) + 2vi ■ SiVi ■ Si 



3vi ■ S1V2 ■ Si — 6n ■ vin ■ SiVi ■ Si + 9n ■ V2n ■ SiVi ■ Si + 3n ■ vin ■ S1V2 ■ Si] 
1) m2G^r /. 4mg 



+ G^l'-^^^^ { 1 + 



C 



(1) rn2GN 



nil 



.)(S?-3(Si.n)')-^(Si.n) 



' oQl 



{Si ■ Si - 35i ■ rvSi ■ n) + (a^i^i))\Sf + vi x Si ■ aj?i) + 1-^2 
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where SlIId is given by (12^ and 5* = ^^e^^^S^^, with S^^ given in (111 I) , with the algebra in 
flASP for the 5 pieces. The squared spin contribution to the EOM in the covariant SSC thus 
read 



c^Si 



[OJ 



J^ + cO J^) X Si + (at?i) X Si) X vi + (Si • vi)u; J^ x vi, (A7) 



-TT - V^s-2 "T ^s2 ; /s oi -r V'ii(i) 

where again 

s2 ^^(1) m^GrN 



..W = 3C^^i.^n(Si-n), (A8) 



and 



^^PN = ^3PN + ^3PN l^yj 

'^SPJV = SCg^^^ [vi(vi • Si) + n X vi (n . (vi x Si)) (AlO) 

-n(n • vi)(Si • vi) - vi(n ■ vi)(Si ■ n)] 

with uj'^pj^ given in ( l56l) . We can now compare the expression in ( l52l) with ( ]A7p and ( ]A9p . 
Using 



(Bi + (^2^;^ X vi) X vi) X Si = iufpj, X Si = -3452^^v2(n ■ Si)(n x Si) (All) 
one can show both spin EOM are indeed identical and the equivalence is thus proven. 

APPENDIX B: CHECKING THE EXTREME MASS RATIO LIMIT 

In this appendix we compare our results to the extremal limit. That is, we consider the 
motion of a test particle with mass m and velocity v, moving in the background of a spinning 
BH with mass M and spin S. In this limit the potential is given by (recall Ces^ = 1 for a 
Kerr BH) 

^- = -^(S^ - 3(S ■ n)^) + ^(S^ - 3(n ■ S)^) - ^(S ■ n)^ 

2^^^v2(S2-3(S-n)2), (Bl) 



4Mr3 
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where we neglect term of order 0{m?) (including the acceleration-dependent part of our 
potential). 

Now we would like to compare this to the effective action generated via 



S^-,nlM,l,/-^'^ (B2) 



by using the harmonic gauge metric 29 1 



ds^ = goodr + 2go^dtd(j) + dli + c?/^2 (B3) 

where the relevant pieces are (with signature (—1,1,1,1)) 

2GnM 2G%M^ GmS^/M - AG%M^ + 3GnS^/M cos(20) 



5'oo = -1 
4 



2G%S^ - 2G%M^ + 4G^^2 cos(2^) 



^4 



2SGNsm^e 2SMG%smH . . 



90^ = + f + 0(G^, S') (B4) 

/y' ly ^ 



„2 GnS^ . 2nf ^. o /omN G7vM(cOs(2^) - 3) ' . , ., 

dll = ^^sm^^ ( -(1 + 3cos(2^)) + "^ ^ ^ ' '- ] d<P' (B5) 



dll = ((r + GnM^ sin^ 9 + 0{GNS))d(P^. (B6) 

For simplicity of comparison, and with no loss of generality, we consider the motion in a 
circular orbit. Notice that we do not expect our results to agree with those in the harmonic 
gauge to all orders in G^, since we are working in a "background field harmonic" gauge. 
However, one can show that the coordinate transformation to the harmonic gauge will have 
no effect to the order we are working in this paper^. Indeed plugging in flB4p into flB2p (after 
changing to our signature convention (1,-1,— 1,-1)) and expanding leads to the result in 



^ Going to the harmonic gauge entails adding a new diagram similar to Fig. [T] with a three graviton 
interaction that follows from Sgf \M- However, one can show that this extra diagram vanishes and our 
result for the S'^ and spin-spin potentials agree with the harmonic gauge to 3PN. That is not the case for 
the spinless part of the potential already at 2PN. We thank Andreas Ross for pointing this out to us. 
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